Quantum magneto-oscillations in a two-dimensional Fermi liquid 



o 
o 

(N 
X> 

IX, 
co 



ed 
-fi 



ed 



I 

-a 
e 

o 
o 



> 

o 
ri 
o 
co 
o 

-I— > 
cd 



i 

fi 

O 
o 



Gregory W. Martin, Dmitrii L. Maslov, and Michael Yu. Reizer* 
Department of Physics, University of Florida, P.O. Box 118440, Gainesville, FL 32611 

(February 2, 2008) 



Quantum magneto-oscillations provide a powerfull tool for quantifying Fermi-liquid parameters of 
metals. In particular, the quasiparticle effective mass and spin susceptibility are extracted from the 
experiment using the Lifshitz-Kosevich formula, derived under the assumption that the properties 
of the system in a non-zero magnetic field are determined uniquely by the zero-field Fermi-liquid 
state. This assumption is valid in 3D but, generally speaking, erroneous in 2D where the Lifshitz- 
Kosevich formula may be applied only if the oscillations are strongly damped by thermal smearing 
and disorder. In this work, the effects of interactions and disorder on the amplitude of magneto- 
oscillations in 2D are studied. It is found that the effective mass diverges logarithmically with 
decreasing temperature signaling a deviation from the Fermi-liquid behavior. It is also shown that 
the quasiparticle lifetime due to inelastic interactions does not enter the oscillation amplitude, 
although these interactions do renormalize the effective mass. This result provides a generalization 
of the Fowler-Prange theorem formulated originally for the electron-phonon interaction. 



Patterns of quantum magneto-oscillations in thermo- 
dynamic (de Haas-van Alphen effect) and transport 
(Shubnikov-de Haas effect) quantities encode three im- 
portant parameters of a Fermi-liquid metal. The pe- 
riod of the oscillations gives the area of the extremal 
cross-section of the Fermi surface, the slope of the tem- 
perature dependence of the oscillations amplitude pro- 
vides the quasiparticle effective mass, and the phase 
shift between oscillations of spin-up and-down electrons 
yields the (renormalized) spin susceptibility. Magneto- 
oscillations studies of the Fermi-liquid (FL) state in 
two-dimensions (2D) date back to the early 70s, when 
semiconductor heterostructures first became available [1] . 
Another surge of the activity in this field, which occurred 
in mid 90s, was stimulated by the discovery of the metal- 
lic state at v = 1/2 [2]. Recently, Shubnikov-de Haas 
oscillations have been used to determine the parameters 
of the "anomalous" metallic state in Si MOSFETs and 
other semiconductor heterostructures exhibiting an ap- 
parent metal-insulator transition in zero magnetic field 
[3-6]. Despite the long and successful history of quan- 
tifying FLs in 3D via magneto-oscillations, this method 
remains controversial in 2D. The primary goal of our pa- 
per is to resolve some of the open issues. 

The first controversy is related to the applicability 
of the current theory of magneto-oscillations to the 2D 
case. The analysis of the experimental data in 2D is 
often based on the premise that the classic result for 
magneto-oscillations in a Fermi liquid for the 3D case, 
known as the "Lifshitz-Kosevich (LK) formula" [7-9], is 
transferrable to 2D upon a trivial change in the electron 
spectrum. The cruciai features of the LK formula, Le., 
its validity for arbitrarily strong interactions (without de- 
stroying the FL) and the fact that the FL parameters 
cntcring the formula are taken at zero magnetic field, 
survive on this premise. The deviations of the observed 
oscillation pattern in stronger fields from that predicted 



by the LK formula are ascribed to oscillations in the ef- 
fective <7-factor [1] and the effective mass [10]. On the 
other hand, there have been warnings that the 3D LK 
formula is non-transferrable to 2D [11,12] for any field 
strength. Hcnce the situation needs to be clarified. The 
second controversy-not specific to 2D-is related to the 
effect of quasiparticle damping. It is often mentioned in 
the literature that any scattering of quasiparticles, elas- 
tic and inelastic, contributes to the smearing of magneto- 
oscillations via the effective Dingle temperature (scatter- 
ing rate) for a given process. Alternativcly, Fowler and 
Prange [13] showed that the electron-phonon scattering 
rate does not appear in the oscillations amplitude due to 
the cancellation of two T— dependent parts of the Mat- 
subara self-energy (cf. also [14]). To the best of our 
knowledge, this cancellation has never been discussed for 
other interactions, including the electron-electron one, 
which is of a primary importance for 2D electron systems. 
The last issue to be addressed in this paper (and not dis- 
cussed previously in the literature) is the effect of inter- 
ference between electron-impurity and electron-electron 
scattering on magneto-oscillations, neglected in the LK 
formula. The theory of interference effects in the ballis- 
tic regime [15], when Tr 3> 1, where r is the electron- 
impurity scattering time (we set h = ks — 1 throughout 
the paper), offers a plausible explanation of the metal- 
lic temperature dependence in the metallic phase of the 
2D metal-insulator transition. Unusual (within the LK 
framework) temperature dependences of the oscillation 
amplitude are also commonly observed in Si MOSFETs 
[6,16], but the proper theory is currently lacking. 

Our answers to these open questions are as follows. i) 
Although it is true that the LK formula does not work 
in 2D at T = and in the absence of disorder, it is 
stili applicable to the situation when finite temperature 
and/or disorder cause the oscillations to be exponentially 
small. ii) The cancellation of the scattering rate term in 
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the Matsubara self-energy is pertincnt to any inelastic 
interaction, including the electron-electron one. Due to 
this cancellation, the scattering rate of inelastic processes 
does not enter the oscillation amplitude. iii) Intcrference 
between electron-impurity and electron-electron interac- 
tions gives a new, TlnT, dependence of the amplitudc's 
argumcnt, that can be interpreted equivalently either as 
a "T— dependent" effective mass or Dingle temperature. 
The functional form of this dependence is the same in 
the diffusive (TV <C 1) and ballistic (TV 3> 1) regimes. 

We limit our analysis to the de Haas-van Alphen effect 
and assume the chemical potential is fixed. The main 
features of the results for the de Haas-van Alphen effect, 
in particular, the T-dependence of the oscillation ampli- 
tude, are commonly expected to apply to the Shubnikov- 
de Haas effect as well, although a rigorous proof of that 
is currently lacking. Assuming a fixed chemical potential 
is not essential for the case of small oscillations, which is 
the focus of this paper (see below). We begin with a brief 
rcminder of how the LK formula is derived in the Lut- 
tinger formalism [8,9]. The key issue here is whether the 
zero-field FL parameters determine uniquely the oscilla- 
tion pattern in a finite (and not small) field. The (Mat- 
subara) self-energy (that encodes ali FL parameters) has 
consists of parts S = S + S osc , where S may contain a 
monotonie (nonoscillatory) dependence on magnetic field 
and S osc oscillates with the field. For electron-electron 
interactions in 3D, |E OSC |/|E | ~ iV~ 3 / 2 , where N is the 
number of occupied Landau levels, whereas the leading 
term in the oscillatory part of the thermodynamic po- 
tential fi falls off as A~ 5 / 2 for S osc = 0. Expanding fi 
in Sosc up to the second order (the first order term van- 
ishes due to the property SCI/5T, = 0), one finds that the 
oscillatory part of S can always be neglected in the semi- 
classical regime (N ^> 1). With this simplification and 
for a momentum-independent self-energy, the amplitude 
of the k th harmonic in Sì is given by 



v = £ cxp (_27rfe[ e , + iMfen,r)r ) j (1) 



U) c 



where e n — 7r(2n + 1)T and uj c = eB/mc. Notice 
that iS is real and does not contain a Constant term. 
The second argumcnt of S emphasizes the fact that 
the temperature enters So in two ways: via the Mat- 
subara frequency and via the thermal distribution of 
electrons and other degrees of freedom. For a generic 
Fermi liquid and in the presence of short-range impuri- 
ties, iS (ie n , T) = ae n + sgne„/2r, so that the effective 
mass is defined as m* = m(l + a). The amplitude then 
assumes a familiar form 



A k = 



2ir 2 kT/uj c 
sinh (2n 2 kT/LU*) 



cxp - 



2^ 2 Tb 



(2) 



where lo* = eB/m*c and Td = l/27rr is the Dingle 
temperature. Momcntum-depcndence of So of the form 
[ìvp (p — pp) results in a change of the effective mass in 



(2) to m* = m (1 + a) (1 + /3) _1 and in multiplying (2) 
by Z s m* /to, where Z s is the renormalization factor. 

In arbitrary dimensionality D, the estimates for the 
ratio of oscillatory to monotonic-in-field parts of the self- 
energy and for the leading oscillatory term in fi change 
to N~ D / 2 and AT-( I3 + 2 )/ 2 , respectively. For D = 2, the 
oscillations in the self-energy are as important as in the 
thermodynamic potential itself [12]. The Luttinger ex- 
pansion at T = Td = breaks down and the LK formula 
is not, generally speaking, valid [12]. The physical rea- 
son is that the ground states of an interacting system 
at B = and in a finite field are not adiabatically con- 
nected in 2D. This fact has been emphasized by recent 
findings that the ground state of a 2D electron liquid is 
not a Fermi liquid even for N ^> 1, but rather a charge- 
ordered state [17]. Nevertheless, an absense of the full 
LK formula in 2D does not preclude a canonical analysis 
of magneto-oscillations, if under more restrictive condi- 
tions, as the FL-behavior is restored at higher energies. 

The power-counting argument for (against) neglccting 
the oscillatory part of the self-energy in 3D (2D) is valid 
at T = T D = 0. If the real and/or Dingle temperatures 
are sufficiently high, i.e., 



2tt 2 (T/u*+T D /u c )>l, 



(3) 



the amplitudes of ali oscillatory quantities, including the 
self-energy, are exponentially small. Neglecting the oscil- 
latory part of the self-energy, the amplitude of the first 
harmonic takes the form 

Ai = (4tt 2 T/lo c ) exp (-2tt [ttT + iS {inT, T)} /u c ) , (4) 

where, due to condition (3), we limited the Matsubara 
sum by the first term eo = ttT. The oscillatory part of 
S results in a correction to Ai which is itself of order 
Ai (with exponential accuracy) . The net contribution to 
fi is of order A 2 , which is of the same order as A 2 for 
S osc = 0. Thus, harmonics with k > 2 are affected by 
the oscillations in S and a 2D analog of the LK formula, 
which includes the sum over ali k, can only be derived in 
a perturbation theory for a specific interaction but not 
for a generic FL. However, the k = 1 harmonic does not 
include S osc and, as long as (3) is satisfied, the analysis 
can proceed as in the 3D case. In what follows, we as- 
sume that (3) is satisfied and the amplitude of the first 
(and only important) harmonic is given by (4). 

Next we discuss whether the quasiparticle relaxation 
rate affeets the amplitude of magneto-oscillations. We 
set Td = temporarily. Suppose that a quasiparticle re- 
laxation rate is measured in a clean Fermi liquid, e.g., via 
electron heat conductivity, with a result that l/r c _ c oc T 2 . 
It seems naturai to assume that the same rate contributes 
also to the Dingle temperature of magneto-oscillations. 
That this is not the case was shown for the electron- 
phonon interaction by Fowler and Prange [13]. Here we 
generalize their arguments for the electron-electron in- 
teraction in 3D and 2D, and then give a general theorem 
for an arbitrary interaction. For a generic FL in 3D, the 
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Matsubara self-energy, up to the quadratic in e n and T 
terms, can be written as 

iS (ie n ,T) = ae n + i(3v F (p - Pf) + 7 i^T) 2 - e 2 n . (5) 

In addition to a direct calculation, the validity of the 
quadratic term in (5) is readily establishcd by noticing 
that upon analytic continuation ie n — > e + iO this term 
gives the correct form for the imaginary part of the on- 
shell self-energy: -ImSj ex (nT) 2 + e 2 [18]. The am- 
plitude of the first harmonic (4) contains iE (inT,T), in 
which the quadratic term vanishes identically. The T 2 - 
terms from higher Matsubara frequencies (legitimately 
considered within this scheme in 3D) increase the ampli- 
tude and cannot be interpreted as "damping." 

In 2D, the integrai over momentum transfers diverges 
logarithmically at the lower limit, changing the behavior 
of Im to E 2 In E, where E = max{e, T}. This change 
does not alter the principal result. Consider the simplest 
case of a contact interaction. To the second order in this 
interaction, the quadratic term in (5) is replaced by 



it (ie n ,T) 



U 2 m T 



co m \n(e F /co m ) . (6) 



Although the sum in (6) does not have an analytic solu- 
tion, it obviously vanishes for e n = . 

To analyze a general case of a finite-range and dynamic 
interaction, including the scrccncd Coulomb one, it is 
convenient to find the imaginary part of the retarded 
self-energy first and then continue back to Matsubara 
frequencies. On the mass-shell, Im T, R is given by 



Im£*(e) = - 



7T 



2(2tt) 



dujF (co) 



coth 



CO 

2T 



tanh 



co 



2T 



where F (co) = J d D qS (co — vf • q) Im V R (co, q) and 
V R (co, q) is the retarded interaction potential. As a func- 
tion of a complex variable z, f(z) = ImE fl (z) has the 
following properties in the upper half-planc: i) ali lines 
Imz = n (2n + \)T are branch cuts on which Re/ is 
continuous but Im/ changes jumpwise; ii) due to the 
fact that tanh(x — ìn(n + 1/2)) = cothx, ali points 
z = in (2n + 1) T are zeroes of f(z). Thus, function f(z) 
is analytic in the band < Imz < nT including the 
point z = inT . Analytic continuation from the real axis 
into this band is legitimate and at z = inT it yiclds the 
Matsubara self-energy E (i^o — ìkT, T) , which is equal 
to zero. Zeroes of / (z) at z = in (2n + 1) T with n > 1 
do not lead to vanishing of £o(*£n, T) for n > 1 because 
those zeroes are separated from the real axis by branch 
cuts and thus are not accessible by analytic continuation. 
The 2D case is special only in that the q— integration rc- 
sults in the Inai factor in F(co) which does not change the 
reasoning given above. In particular, for a dynamically 
screened Coulomb interaction in 2D, F (co) oc co In \co\ and 
stili È (ìttT,T) = 0. As this result does not depend on 



the particular form of the interaction, it can be viewed 
as a generalization of the Fowler-Prange theorem. No- 
tice that in 3D the Fowler-Prange theorem is of limited 
applicability because nothing prevents one from consid- 
ering k > 1 and lower values of T + To, when the effect 
of T,o(ie n> o, T) needs to be taken into account. In 2D, 
one is bound to consider only T, (ie n =o,T) within the 
Luttinger approximation. 
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FIG. 1. (a) The interference correction to the self-energy. 
(b) The vertex correction is assigned to either one of the 
vertices in (a) because the self-energy arises as insertions 
into the thermodynamic potential (closed loops). (c) Sin- 
glet/triplet-channel contributions to the effective potential. 

Finally, we discuss the effect of interference between 
electron-electron and electron-impurity scattering on 
magneto-oscillations, extending the analysis of the in- 
terference corrections to the self-energy in 2D from the 
diffusive (Tt « 1) [19] to the ballistic (TV > 1) limit. 
The general form of the interference correction to the 
Matsubara self-energy is (see Fig. 1) 



4 nt (^n,p) =-2T f^oV(ico m ,q) (7) 

e n {u m -e n )>0 J ( 27F > 



(2nY 

xT (ico m , q) G (ie n - ico m , p q) 



where G (ie n , p) = (is n — £ p + zsgn£„/2r) and the ef- 
fective interaction V = V s + 3V t contains the contribu- 
tions from both singlet and triplet channels [15] 

V- 1 (ico m , q) = (2ne 2 /q + F^/v)^ 1 - li (ico m , q) ; 
Va 1 ( iuJ m, q) = vjF^ - LI (ico m , q) , 

where v — m/ir. The factor of two in (7) ac- 
counts for two possibilities of including the vertex cor- 
rection in Fig. la. The Fermi- liquid constants F9 and 
F® [20] determine the rcnormalizcd charge and spin- 
susceptibilities, respcctively. The vertex T(ico m ,q) = 

1) _1 reduces to T = 

-1/2 



co m \ 2 + (qv F ) 2 



(\J (\u) m \ t + l) 2 + (qv F TY 

(Drq 2 + \cOm\Ty 1 and T = r" 1 

in the diffusive and ballistic limits, respectively. The gen- 
eral form of the (small q) polarizarion operator 

LI (ico m ,q) = -v[l- \co m \rT (ico m , q)} (8) 
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reduces to II (iu> m , q) = —vDq 2 / (Dq 2 + \uj m \) in the dif- 
fusive limit, whcrc D = v f t/2, and to II (iu) m ,q) = 

— v ^1 — |w TO | I \j (v F q) 2 + uj^j in the ballistic onc. 

Omitting the details of lengthy but straightforward cal- 
culations, we give just the result for the self-energy valid 
to logarithmic accuracy 

? 4 nt (inT, T) = -Tln (e F /T) Q (Tt, F°) /2s f t, 

where Q (Tr, F°) = g p (Tt) + [3F°/ (l + F°)] g a (Tt) 
and g p / a (x) are slowly varying functions which inter- 
polate between the diffusive and ballistic regimes. The 
limiting valucs of g p / a ( x ) are &s follows: g p (0) = 1, 
g p (x > 1) = 3/2, g„ (0) = 1, g a (x > 1) = 1/2. Apart 
from the numcrical cocfficients, the T— dependencc of 
iSg nt (7rT, T) is the sanie in the diffusive and ballis- 
tic regimes. In that sense, the behavior of the self- 
energy is sumilar to that of the tunneling density of 
states [21]. Notice that the interference correction to 
the scattering rate in the ballistic regime | Im | ~ 
(T / e j?t) In (e f /T) is smaller than the scattering rate 
in a clcan FL, |lmE^| int - {T 2 / e F )\n(e F /T), in the 

parameter (Tr)~ <C 1. However, due to the can- 
ccllation of the T 2 lnT term in iT,o(iirT, T) discussed 
carlicr in this paper, the interference correction is the 
main nonlinear T— dependent term in the Matsubara self- 
energy, leading to a modification of the LK-formula. The 
T In T— dependence of the self-energy can be interpreted 
as a logarithmic T— dependence of the effective mass. 
Following this interpretation, the interference effect leads 
to a replacement of the effective mass in the argument of 
the exponential in Ai by 

Notice that the effective mass is reduced by the singlet- 
channel interaction but enhanced by the ferromag- 
netic (F® < 0) interaction in the triplet channel. The 
In T— dependence of the effective mass is a characteristic 
feature of the "marginai Fermi liquid" model [22,23]. 

Equivalcntly, the nonlinear T— dependence of 
iT,o (inT, T) may be interpreted as a T— dependent Din- 
gle temperature: 

T D (T) =T D (1- (T I e F ) In (e f /T) Q (Tt, F )) . (10) 

One of the empirical procedures used in [6] to account for 
the observed deviations from the LK formula was to as- 
sume that the effective Dingle temperature has the same 
T— dependence as the zero-field resistivity. Comparing 
(10) with the result for the interference correction to the 
resistivity [15], we see that, although this recipe is not 
precise, it has some theoretical justification: the general 
structure of the Tp, (T) and p (T) depedences is similar 
in the ballistic regime, except for a factor of In (e f /T) 
present in Tp (T) but not in p (T). 
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